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Abstract. We show that the relic abundance of the minority component of asymmetric dark 
matter can be very sensitive to the expansion rate of the Universe and the temperature of 
transition between a non-standard pre-Big Bang Nucleosynthesis cosmological phase and the 

y—i standard radiation dominated phase, if chemical decoupling happens before this transition. 

In particular, because the annihilation cross section of asymmetric dark matter is typically 

r^ larger than that of symmetric dark matter in the standard cosmology, the decrease in relic 

^sO density of the minority component in non-standard cosmologies with respect to the majority 

component may be compensated by the increase in annihilation cross section, so that the 
annihilation rate at present of asymmetric dark matter, contrary to general belief, could be 

-^ larger than that of symmetric dark matter in the standard cosmology. 
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1 Introduction 

The nature of dark matter (DM) is one of the fundamental problems of physics and cosmol- 
ogy. Particles with weakly interacting cross sections and masses in the few GeV to 10 TeV 
range, WIMPs (Weakly Interacting Massive Particles), are among the best motivated DM 
candidates. DM particle candidates, such as WIMPs (but also sterile neutrinos and axions) 
are produced before Big Bang Nucleosynthesis (BBN), an epoch from which we have no data. 
BBN is the earliest episode (finishing 200 seconds after the Bang, when the temperature of 
the Universe is T ~ 0.8 MeV) from which we have a trace, the abundance of light elements 
D, 4 He and 7 Li. In order for BBN and all the subsequent history of the Universe to proceed 
as usual, it is enough that the earliest and highest temperature during the last radiation 
dominated period, the so called reheating temperature Tru, is larger than 3.2 MeV [1]. 

The argument showing that WIMPs are good DM candidates, many times called the 
"WIMP miracle" , is more than 30 years old [2] . The density per comoving volume of non 
relativistic particles in thermal equilibrium in the early Universe decreases exponentially with 
decreasing temperature, due to the Boltzmann factor, until the reactions which change the 
particle number become ineffective. At this point, when the annihilation rate becomes smaller 
than the Hubble expansion rate, the WIMP number per comoving volume becomes constant. 
This moment of chemical decoupling or freeze-out happens later for larger annihilation cross 
sections a, which produces smaller WIMP densities. If the Universe is radiation dominated 
during decoupling and there is no subsequent change of entropy in matter plus radiation, the 
present relic density is 



which for weak order a ~ G^m?, gives the right order of magnitude of the DM density 
(and a temperature TV . ~ m x /20 at freeze-out for a WIMP x of mass m x ). 

The "WIMP miracle" argument and the standard computation of relic densities rely on 
assuming that radiation domination began before the main epoch of production of the relics, 
that the entropy of matter and radiation has been conserved during and after this epoch, 
that WIMPs are produced thermally, i.e. via interactions with the particles in the plasma, 
and that there is no significant asymmetry between the WIMP particles and antiparticles. 
With these assumptions chemical decoupling happens at Tj ~ m x /20, thus all WIMPs 
with m x > 80 MeV decouple at temperatures higher than 4 MeV, when the content and 
expansion history of the Universe may differ from the standard assumptions. 

The relic density, see e.g. Ref. [3] (and also the relic velocity distribution, see e.g. 
Ref. [4]) before structure formation of WIMPs (and other DM candidates, e.g. sterile neu- 
trinos [5]) depends on the characteristics of the Universe (expansion rate, composition, etc.) 
before BBN. If these particles are ever found, they would be the first relics from the pre-BBN 
epoch that could be studied. Thus we will want to extract as much information about the 
Universe at the moment these particles decoupled as we can. 

Here we present a new potentially detectable effect that non-standard cosmologies may 
have on asymmetric dark matter. 

2 Asymmetric dark matter 

The idea of asymmetric DM is almost as old as the "WIMP miracle" argument: if DM 
particles and antiparticles have an asymmetry similar to the baryonic asymmetry this could 
explain why the baryonic and DM relic density are similar, i.e. differing by a factor of a 
few and not by many orders of magnitude. In 1985, S. Nussinov [6] pointed out that if the 
asymmetry of technibaryons and usual baryons could be similar in the early Universe, their 
present number relic density would also be similar, which would mean that the ratio of their 
relic densities would be given by the ratio of the DM and baryon masses. It was then largely 
assumed then that the DM relic density was the critical density, thus SIdm/^Ib — 100 (for 
£Ib — 0.1) would imply a ratio of the lightest neutral "Technibaryon (TB)" and nucleon 
mass to be rriTB/GeV ~ 100 which was phenomenologically acceptable at the time. This 
argument does not exactly hold now, because the ratio of dark and visible matter relic 
densities is only about a factor of 5, but different versions of Technicolour models have been 
a fertile framework for 100 GeV-TeV mass asymmetric DM candidates (see e.g. [7]). 

The first model unrelated to Technicolour for "light" asymmetric DM was proposed 
in 1986. Gelmini, Hall and Lin [9] produced several models for 5 to 10 GeV mass DM 
candidates, called "cosmions" (X), which could influence the physics of the Sun. In one 
of them the baryon minus "cosmion" number {B — C) is assumed to be conserved and a 
common origin of an asymmetry in both B and C numbers, AB = AC insures that the 
baryon and DM asymmetries are identical, so that Qx/^B = mx/GeV~ 5 to 10 (which was 
then considered enough just to account for the "galactic dark matter"). The same basic idea, 
i.e. the conservation of a linear combination of baryon and DM particle number and a shared 
origin of the baryon and dark matter asymmetry was subsequently realized in many models 
for heavier candidates, in the 10's to 100's GeV mass range (see e.g. Ref. [8]), until Kaplan, 



Luty and Zurek in 2009 [10] applied it again to candidates in the 5 to 15 GeV mass range, 
in a paper which generated renewed interest in asymmetric DM models (for a recent review 
see e.g. [11]). Light asymmetric DM particles not unlike "cosmions" were recently proposed 
as a means to solve the current discrepancy in the predicted composition of the Sun between 
helioseismological data and the revised Standard Solar Model [12]. 

The ideas we develop in the following apply to asymmetric DM WIMPs of any mass. 
We are not concerned here with the origin of the asymmetry. We assume only that one was 
generated before the chemical decoupling of the DM particles x an d antiparticles x, so that 

Y X -Y X = A, (2.1) 

where Y x = n x /s, Y x = n x /s, n x and n x are the respective relic number densities, s is 
the entropy density, s = (2ir 2 / 45) <?*T 3 , dominated by the relativistic degrees of freedom g* 
and A is a constant that characterizes the asymmetry. Here we take A positive so x an d X 
are respectively the majority and minority components of the DM at present. (Note that 
our asymmetry parameter A is called 2a in Ref. [13] and C in Ref. [14]). Moreover, in the 
following we will assume that x an d X account for the whole of the DM, i.e. 

^x + ^x = ^dm, (2.2) 

although the arguments can be easily changed if they account for only a fraction of the DM. 
The evolution of the equilibrium values of Y x , Y x as function of x = rn x /T are shown 

in Fig. l.a. The equilibrium number densities n x and n x in presence of an asymmetry 
differ by the chemical potential fi x (in equilibrium fi x = —fi x ) 
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Here, m x is the mass of x an d X, and g x is the number of internal degrees of freedom of x 
and of x separately. For simplicity in the following we take g x = 1 (which assumes x an d X 
are conjugate complex scalar fields). A different choice of g x (e.g. g x = 2 if x and X are two 
conjugate Dirac fermions) would not affect the results significantly. The chemical potential 
can be written in terms of the asymmetry A by substituting the equilibrium number densities 
into Eq. (2.1), 

3 

nf2( % = 0) (W T - e~^ T ) = g x {^) ' e~ m - /T (e^ T - e~^ T ) = As. (2.5) 
and solving this quadratic equation for e^ x ' T , to get 



In the presence of the asymmetry A, the annihilation of the majority component is 
considerably reduced after Y x reaches the value A at xa, Y x (xa) — A. For x > xa the 
equilibrium number density per comoving volume of the majority component becomes almost 
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Figure 1. l.a (left) Evolution of the equilibrium abundances Y E ® — n E< ^ / s for the majority, %, and 
minority, %, dark matter components, as function of x — rn x /T for A = Y x — Y x = 4.05 x 10~ 12 and 
m x = 100 GeV. Also shown is the equilibrium abundance of symmetric dark matter (A=0). l.b (right) 
Equilibrium annihilation rates of x an d X, ^ x and ^ x respectively, for three increasing values of 
the annihilation cross section (cr xx v), 9.5 x 10~ 9 GeV~ 2 , 9.0 x 10~ 7 GeV" 2 and 5.0 x 10~ 6 GeV~ 2 , for 
the lower (black), middle (blue) and higher (red) lines respectively. 



constant while that of the minority component decreases faster than in the symmetric A = 
case as x increases. It is easy to see why this is so by considering the annihilation rate of x 
and x, T x and T<> respectively, 
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in equilibrium, whose evolution as function of x is shown in Fig. l.b. Here (a xx v) is the 
thermally averaged XX annihilation cross section and we are assuming that these annihilations 
are the only processes that can change the number of these particles. 

For x > xa the x an d X interaction rate in equilibrium become respectively larger and 
smaller than in the symmetric A = case because n x is much larger and n x is smaller than 
in the ^4 = case. 

The annihilation of each DM component x an d X ceases when their respective interaction 
rates become smaller than the expansion rate of the Universe H, which happens at their 
respective decoupling or freeze-out, Xf for x & n d x~f for %, defined by 



T x( x fo) = H(x fo ), T x (x fo ) = H(x fo ). 



(2i 



It is clear from Fig. l.b that x decouples earlier than x, i.e. Xf < Xf . At their respective 
freeze-out, each component acquires its relic density, fixed in comoving volume. However, 
little changes in the abundance of the majority component after it freezes-out, since already 
at xa < Xfo the annihilations had become very suppressed. In the standard cosmology the 
relic density of the minority component, which was computed for the first time in Ref. [13] 
and recently in Ref. [14], is exponentially small with respect to the majority component 
density. This is why there is no asymmetric DM annihilation at present. 
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Figure 2. Expansion rate of the Universe H for the standard cosmology, STD (solid red line) and 
for the kination, K (dashed blue line) and scalar-tensor, ST (dot-dashed green line) cosmologies with 
a transition temperature T tr = 5 MeV to the standard cosmology. 



Here we will show that if the chemical decoupling of the minority component hap- 
pens during a non-standard cosmological phase, the relic abundance of this component of 
asymmetric dark matter can be much larger than in the standard cosmology. Thus, if the 
annihilation cross section of the dark matter candidate is known, the annihilation rate at 
present, if detectable, could be used to test the Universe before Big Bang Nucleosynthesis, 
an epoch from which we do not yet have any data. In particular, because the annihilation 
cross section of asymmetric dark matter is typically larger than that of symmetric dark mat- 
ter in the standard cosmology, the decrease in relic density of the minority component with 
respect to the majority component may be small enough that the annihilation rate at present 
of asymmetric dark matter, contrary to general belief, could be larger than that of symmetric 
dark matter in the standard cosmology. 



3 Non-standard pre-BBN cosmologies 

Many viable non-standard pre-BBN cosmological models have been proposed (see e.g. [15] 
for a review). Usually non-standard cosmological scenarios contain additional parameters 
that can be adjusted to modify the WIMP relic density. However these are due to physics 
that does not manifest itself in accelerator or DM detection experiments. 

In the standard (STD) cosmology the Universe is radiation dominated, thus the expan- 
sion rate of the Universe H depends on the temperature of the radiation bath T as H ~ T 2 , 
but in non-standard models H can decrease slower, H ~ T 12 , or faster, H ~ T 3 as T de- 
creases. This is what happens in "scalar-tensor" [16, 19] and "kination" [20] cosmological 
models, respectively. In both of these models the entropy of matter plus radiation is con- 
served (and it is dominated by the radiation component), thus the usual relation between T 
and the scale factor of the Universe a, namely a ~ T , is the same as for a radiation domi- 
nated Universe. In models in which the potential energy of a scalar field oscillating around 
its true minimum while decaying is the dominant component of the Universe just before 



BBN, the expansion rate decreases even faster, H ~ T but there is entropy creation, thus 
T ~ a~ 3 ' 8 . These models produce the largest departure of symmetric WIMP relic density 
from the standard density (see e.g. [3]) and we could expect their effect on asymmetric DM 
to be very large too. However here for simplicity we concentrate on models in which there is 
no entropy creation. 

Scalar-tensor (ST) theories of gravity [16, 19] incorporate a scalar field coupled only 
through the metric tensor to the matter fields. In many of these models the expansion of 
the Universe drives the scalar field towards a state where the theory is indistinguishable 
from General Relativity at a transition temperature Tt r , but the effect of the scalar field 
changes the expansion rate of the Universe at earlier times, either increasing or decreasing 
it. Theories with a single matter sector typically predict an enhancement of H before BBN. 
In Ref. [19], for large temperatures T > 2j r the H is enhanced with respect to the standard 
expansion rate #std by a factor fy ~ 2.19 x 10 14 (To/T) 0,82 (To is the present temperature 
of the Universe), thus H$t ~ T ■ At T tr , f^ drops sharply to values close to 1 at T < T tr 
before BBN sets is. 

Kination [20] is a period in which the kinetic energy p^ ~ 2 /2 of a scalar field <fi 
(maybe a quintessence field) dominates over the potential energy so Ptotai — <f/2 ~ o~ 6 . 
Since T ~ a -1 as usual, H ~ y'ptotai ~ T 3 . The contribution of the 4> kinetic energy 
to the total density is usually quantified through the ratio of <j> -to-photon energy density, 
V4> = P<t>l 'Pi at T ~ 1 MeV so that at higher temperatures H ~ ^/ly ) (T/lMeV)HsTD- Notice 
that at T ~ 1 MeV, i.e. during BBN, the quintessence field cannot be dominant, thus rjj, < 1. 

Fig. 2 shows the standard (STD) expansion rate, and the scalar-tensor (ST) and kina- 
tion (K) expansion rates for a transition temperature (Tt r is the temperature at which the 
cosmology becomes standard) Tt r = 5 MeV. 

4 Minority component relic density calculation 

The relic densities of DM particles x an d anti-particle x are calculated by solving their 
respective Boltzmann equation (for asymmetric DM, x an d X are n °t self-conjugate, x 7^ X)- 
We assume that close to the moment of decoupling, the only reactions that change the number 
of DM particles and antiparticles are annihilations of XX pairs into Standard Model particles 
and the inverse process of pair creation. With this assumption we have 

^ + 3Hn x = ^ + mn t = -{a xt v){n x n t - n^nf). (4.1) 

In terms of the dimensionless quantities Y x , Y x and x, and treating g* as a constant, the 
Boltzmann equations become 

dY x _ dY x _ -(a xx v) 2ir 2 m\ y eq y eq, (d<? , 

dx ~ dx ~ H 45 5 *^ 1 x * * * h [ ' 

where (see Eqs. (2.3) and (2.4)) 
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In the following we take g+ = 90 for simplicity, a good approximate value for decoupling 
temperatures above the QCD phase transition. Eq. (4.2) implies d(Y x — Y x )/dx = 0, so 



the asymmetry A defined in Eq. (2.1) is a constant. It is clear that A is constant because 
we are only considering interactions of the form XX *~> P' s where p's are Standard Model 
particles, which leave the difference between the co- moving number density of the particle 
and anti-particle unchanged. We are assuming that the asymmetry has been produced prior 
to the epoch that we are considering. 

Using Eq. (2.1), we finally obtain the equations we wish to solve 
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4.1 Analytical solutions 

It is useful to obtain approximate analytical solutions to Eq. (4.5) and compare them with 
numerical solutions. Fig. l.b makes it clear that in period between xa, when Y x becomes 
practically constant and the freeze-out Xf of the minority component, Y x decreases expo- 
nentially. Until x freezes out, the densities of both DM components are those of equilibrium 
(for x it is very close) and after, for x > xt , the production term in the Boltzmann equations 
is suppressed compared to the annihilation term (because Y x <C Y x ) and so we ignore it. 
Then, for x > xj , Eq. (4.5) becomes 
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Integrating this equation from xj to oo (a good approximation for the value of x at present) 
we get 
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For asymmetric DM the ratio Y x /Y x increases immediately after x decouples, as Y x 
decreases. The decoupling of x is n °t instantaneous, as assumed when defining xf , and 
annihilations still continue for a while after Y x departs from its equilibrium value. This can 
be seen e.g. in Figs. 9 and 10 where the solid black lines showing the numerical solutions for 
Y x decrease after the freeze-out of the minority component at Xf > xa- So we keep only the 
x — > oo term on the left hand side of Eq. (4.7) and obtain 
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Using Eq. (2.1), i.e. Y x (x — > oo) = A + Y x (x — )■ oo) we get for the majority component 



Y x (x -»■ oo) = A 
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4.1.1 Standard cosmology 

In the standard cosmology, the expansion rate of the Universe is 
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(4.10) 



where Mp is the reduced Planck mass, and g* is the effective number of relativistic degrees 
of freedom. Expanding the annihilation cross section in powers of the relative velocity, v, 
and then taking the thermal average, the annihilation cross section can be written in term 
of constants a and b as 

(a xx v) = a + 3bx~ l + 0{x~ 2 ). (4.11) 

where the a term is dominant for s-wave XX annihilation, and a = for p-wave annihilation. 
Using Eqs. (4.10) and (4.11), we solve the integral in Eq. (4.8) with the result 



Y? TD (x -> oo) = A \ exp 
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where A is defined as A = — [(ds/dx)/3H] x=1 
and mioo = m x /100 GeV. From Eq. (4.9), we get 
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4.1.2 Kination 



At the transition temperature T& at which the Universe becomes radiation dominated after 
the kination phase, the expansion rate of the Universe during kination, Hx, coincides with 
the standard one, HxiTtr) — HsTD{Ttr) which allows us to fix the constant multiplying T 3 
in Hk and find 



H K {T) 
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(4.14) 



MpT tr V 90" 

This corresponds to having the kination parameter ^/ff^ = lMeV/Tj r , so for T% T > 1 MeV, 
Tij, < 1 during BBN. 

Using Eq. (4.14) in Eq. (4.8) we calculate the relic abundance of x due to kination, 
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(4-15) 

where x% T = m x /Tt r , and A is defined as before. An expression for the relic abundance of % 
can be found by using Eq. (2.1), i.e. Y x (x — > oo) = A + Y x (x — > oo). 

4.1.3 Scalar-tensor model 

As mentioned above, in this model [19] the expansion rate of the Universe is given by 

Hst(T) = U(T)H S td(T) (4.16) 

where /</>(T) ~ 2.19 x 10 14 (Tq/T) ' is a temperature dependent factor which transitions 
very fast to a value of 1 at the transition temperature Tt r , and Tq is the present temperature 



of the Universe. In terms of x, this factor is f^,(x) — 9.65 x 10 3 (GeV/m) ' x 082 . Thus, for 
transition temperatures from 1 MeV to a few GeV, the Hubble parameter decreases suddenly 
by a few orders of magnitude at Tf r , when the change to the standard cosmology happens. 
Due to this sudden drop in the H, \ an d X start re-annihilating and they freeze-out again 
shortly after. Usually % and \ freeze-out in the re-annihilation phase before they can reach 
equilibrium and, in this case, we can still neglect the production term in the Boltzmann 
equations (again because Y x S> Y x y ), and Eqs. (4.8) and (4.9) still hold. Using Eq. (4.16) 
in Eq. (4.8), we calculate the relic abundance of X i n the scalar-tensor cosmology, 



Y? T (x -)• oo) = A I exp 



(4-17) 

where xt r = m x /Tt r , and A is defined as before. We have taken 6 = for simplicity. As 
before, an expression for the relic abundance of x can De found by using Eq. (2.1), i.e. 
Y x (x — > oo) = A + Y x (x —f oo). 

AAA Estimation of the X freeze-out temperature 

To find Xf we use the freeze-out condition Eq. (2.8) for %, using in each case the correspond- 
ing H written above and the reaction rate of x, T^ considering that until freeze-out n x and 
n x closely track their equilibrium values, n x and n x given in Eqs. (2.3) and (2.4), so we 
can approximate T x by 

I* ~ (a xx v)n^^ = 0K*/ T (4.18) 

Using Eqs. (2.5) and (2.6) we get an expression for T x as function of the asymmetry A. Using 
this form of T x , the X freeze-out condition Eq. (2.8) becomes a transcendental equation for 
Xf , which we can solve numerically. After solving for Xf , we use the expressions calculated 
above to estimate the relic abundance of x an(1 X- 

4.2 Numerical solutions 

The Boltzmann equations (4.4) and (4.5) can be solved numerically for each of the cosmolo- 
gies we consider. Here we show examples of the evolution of Y x and Y x as a function of x in 
each cosmology and later compare the numerically calculated relic density to the analytical 
results. In all of the following we present results for pure s-wave XX annihilation. 

The final relic density of the minority component depends strongly on the relation 
between its freeze-out temperature or Xf and the value xa at which the asymmetry becomes 
important and the minority and majority densities Y x and Y x become very different, 
given by Y x (xa) — A (see in Fig. 2. a). 

One can distinguish three cases. The first is when xj < xa, i-e. the decoupling of X, 
and also x since in this case xj = Xf , happens before the split of Y x and Y^ , then 
Y x ~ Y x at freeze-out, and so the present relic abundance of x an d X are practically the 
same. This case is very similar to that of symmetric DM. 

The second case is when xj S> xa, i-e. when the decoupling of x occurs well after the 
split between Y x ^ and Y^ , so that Y x (xf ) is exponentially smaller than Y x ~ A, thus 
the x present relic density is completely negligible when compared to the x density. 

The last case, in which the final n x is smaller but not much smaller than n x , happens 
if Xf > xa but Xf ~ xa, i-e. if the decoupling of X occurs shortly after the split between 
n x and n x . Given a particular asymmetry A, the range of cross sections that satisfies this 
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Figure 3. Values of m 100 A for which x an d X constitute all of the DM, as function of the \X 
annihilation cross section, in the standard cosmology (blue contour in both panels) and 3. a (left) 
kination and 3.b (right) scalar-tensor models. The solid gray contour to the farthest right in each plot 
corresponds to a transition temperature of T tr = 4 MeV, and the dashed gray contours progressively 
to the left have higher values of T tr approaching the standard cosmology freeze-out temperature (solid 
blue contour). The contours depends very slightly on m x for the standard cosmology, and if m x > 100 
GeV also for the other two (for m x ~ 10 GeV and T tr close to 4 MeV, the contours change appreciably 
with respect to those shown here- see Fig. 5). 



relationship is narrow, approximately within a factor of three. We choose A such that \ and 
X constitute all of the DM, Eq. (2.2). In the standard cosmology the necessary annihilation 
cross section is very close to that of symmetric DM satisfying the same condition, and as soon 
as the cross section departs by a factor of a few from this value, the relic minority component 
density becomes exponentially small. In the kination and scalar-tensor cosmologies, for a 
fixed transition temperature, the required values of the annihilation cross section in this third 
case can be much larger than the standard symmetric cross section (although in each case 
the values must be in a reatively narrow range, within a factor of a few). This is illustrated 
in Fig. 3, which presents the contours of asymmetry value uiioqA (recall mioo = m x /100 
GeV) as a function of the XX annihilation cross section for which x an d X niake up all of 
the DM (Eq. (2.2) holds). As exemplified in Figs. 4 and 5, the miooA contours are almost 
independent m x in the standard cosmology and if m x > 100 GeV in the other two. The 
contours of Fig. 3 do not apply if m x ~ 10 GeV and T tr close to 4 MeV (see Fig. 5). We use 
&DM = (Px + Px)/Pci with Qdm = 0.11/i~ 2 , p c = 1.05375 x 10~ 5 h 2 cm -3 , p x = m x Y x so and 
p x = m x Y x so where so = 2889.2 cm -3 is the present entropy density of the Universe [22]. 
With these values of the various parameters, the miooA contours satisfy the condition 



mioo(^ x + ^x)= 4 - olxl ° 



-12 



(4.19) 



For the standard cosmology, there is only one contour that satisfies Eq. (4.19) (the left most 
contour, shown in blue). In kination and scalar tensor models the contour satisfying the 
same condition depends on the transition temperature Tj r , which in relevant models can take 
values between somewhat below 4 MeV and the standard cosmology freeze-out temperature. 
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Figure 4. Contours of miQoA for which \ an d X make up all of the DM as a function of the \X 
annihilation cross section for m x = lOGeV (green), m x = 100 GeV (red), and m x — 500 GeV (blue) 
in the standard cosmology. 



Thus, there is a range of contours which satisfy Eq. (4.19) for the kination and scalar-tensor 
models, as shown (dashed gray lines) in Fig. 3. a and Fig. 3.b respectively. Values of A and 



xx l 



above to the left of the contours would lead to a relic density larger than the DM 
density. Given a particular contour, values of (cxx v ) an d A below it lead to a relic density 
of x an d X smaller than the DM density. 

For each DM particle mass m x there is a maximum value of A, corresponding to the 
horizontal part of each contour, which satisfies Eq. (4.19) with Y x + Y x = A, which (since 
Y x — Y x = A) means that Y x = A and Y x is negligible. Thus the horizontal part each contour, 
common to all of them, corresponds to the second case mentioned above, where xt 3> xa- 
Only in the curved and vertical parts of the contours can Y x be non-negligible with respect to 
Y x . As the asymmetry A becomes smaller, the contribution of x to the relic density increases, 
and when A=0 we recover the symmetric case. Eq. (4.19) and the definition of A imply that 
for values of mioo^l. ranging from to 3.3 x 1CP 12 , the relic density of x must be within a 
factor of 10 of the density of x- The curved part of the contour corresponds to the third 
case mentioned above, where Xf is larger but very similar to xa, for which the asymmetry 
is large but the relic density of x is only a few orders of magnitude smaller than the density 

of X- 

As already mentioned the contours in Fig. 3 are only slightly affected by changes in m x 
in the standard cosmology, and if m x > 100 GeV also in kination and scalar-tensor models. 
Fig. 4 shows the almost overlapping m\ooA contours for which x an d X constitute the whole 
to the DM in the standard cosmology for m x of 10 GeV (red), 100 GeV (green), and 500 
GeV (blue). The analytic solutions predict this behavior. For example, one can calculate 
the minimum annihilation cross section for a mass m x , which occurs when mioo^- - > 0. For 
small mioo^) the exponential term in the analytic solutions (Eqs. (4.12) and (4.13)) can 
be expanded as exp[AXa/xf ] ~ 1 + AXa/xf . With this approximation, and substituting 
the analytic solutions into Eq. (4.19), we see that 2miooXf /(\a) = 4.01 x 10 -12 . Recalling 
that A = 3.0 x 10 21 mioo GeV and using x ~ 20, as in the symmetric case, we calculate a 
minimum cross section a ~ 3.3 x 10~ 9 GeV~ , independent of m x . 

Figs. 5. a and 5.b show the same as Fig. 4 but for the kination and scalar-tensor mod- 
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Figure 5. Same as Fig. 4 but for 5. a (left) kination and 5.b. (right) the scalar-tensor model. The 
solid contours correspond to a transition temperature of T tr — 4 MeV and the dashed ones to T tr 
close to the standard cosmology freeze-out temperature. In 5.b the dashed lines are superimposed. 



els, respectively. The solid contours correspond to models with a transition temperature of 
T tr = 4 MeV, and the dashed contours to T tr close to the standard cosmology freeze-out tem- 
perature. The figure shows that for m x = 10 GeV the tuioqA contours change considerably 
with respect of those for m x > 100 GeV only if Tt r is low. 

Notice that the value of miooA in the horizontal part of the contours is unaffected by 
changes in m x in all three models, Figs. 4 and 5. Clearly, this must be so because there 
Y x = A, and Y x is negligible, thus the fixed DM density depends only on the product m x A. 

We mentioned three cases above in which either Xf < xa, Xfo ^> %A or xj is larger but 
very close to xa, so that the minority relic number density n x is either equal to, negligible 
or just a few orders of magnitude smaller than the majority relic density n x , respectively. 
Figs. 6, 7 and 8 demonstrate these three cases in the standard cosmology, kination and 
scalar-tensor models respectively. In all three figures the bottom left panel shows the regime 
we are interested in, where the minority component relic density n x is only a few orders of 
magnitude smaller than the relic density of the majority DM component n x . 

5 Comparison of analytic and numerical relic abundance solutions 

Figs. 9 and 10 show that the analytic solutions we found above agree very well with the 
numerical solutions in all three cosmological models we examined. In the examples presented, 
both solutions agree within a factor of two. In both figures m x =100 GeV and A = 4.01 x 
10~ 12 but the cross section for each particular cosmological model is chosen so that xj is 
not much larger than xa, because we wanted to show examples where the relic density of 
the minority component is only a few orders of magnitude smaller than that of the majority 
component (see the figure captions). 

We have ascertained that the analytical solutions for the relic abundances Y x and Y x 
after \ freeze-out reproduce the numerical solution to within a factor of four while Y x /Y x > 
10~ 10 (smaller ratios are completely irrelevant). As this ratio increases, the agreement is 
better, within a factor of two. 
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Figure 6. In 6. a (top left) the equilibrium annihilation rate over {cr x ^v) (black lines) is compared with 
the expansion rate in the standard cosmology over (a xx v), for three values of {a xx v), 1 x 10~ n GeV -2 
(top, red, line), 1 x 10 -8 GeV~ 2 (middle, blue, line) and 1 x 10~ 5 GeV~ 2 (lower, green, line). 6.b (top 
right), 6.c (bottom left) and 6.d (bottom right) show the evolution of Y x (dashed colored line) and Y x 
(solid colored line) compared to their equilibrium values (dashed and solid black lines) for the three 
values of (<J xx v) just mentioned, respectively. In 6.b Y x ~ Y x , and the red lines are superposed. In 
6.c Y x (dashed blue line) closely follows equilibrium, and Y x (solid blue line) results less than 3 orders 
of magnitude smaller than Y x . In 6.d both \ and x closely follow their equilibrium abundances. 



6 Present asymmetric DM annihilation rates 

Given a particular value of the XX annihilation cross section and almost identical values of 
the asymmetry, very different minority over majority asymmetric DM density ratios may 
result in the three pre-BBN cosmologies we explored, if the chemical decoupling happens 
during these different cosmological phases. Thus, if the annihilation cross section of the DM 
candidate is known, the annihilation rate at present, if detectable, could be used to test 
the Universe before BBN, an epoch from which we do not yet have any data. In particular, 
because the annihilation cross section of asymmetric DM can be larger than that of symmetric 
DM in the standard cosmology, the decrease in the minority component relic density may 
be compensated or even overcompensated by the increase in annihilation cross section so 
that the annihilation rate at present of asymmetric DM, contrary to general belief (see e.g. 
Ref. [12]), could be non- negligible and even larger than that of symmetric DM in the standard 
cosmology. Assuming a "neutralino like" DM candidate in the standard cosmology (for which 
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10- 6 GeV~ 



Same as in Fig. 6 but for kination models and (cr xx v) equal to 1 x 10 GeV , 1 x 
and 1 x l(T 4 GeV~ 2 . 



particle and anti-particle coincide) the standard annihilation rate at present can be written 



as 



■^STD 
sym 



-(a se i f v){p 2 DM /m 



(6.1) 



where the 1/2 factor corresponds to identical particles annihilating, pdm = &dmPc is the 
DM energy density, {a se ifv) is the velocity averaged self-annihilation cross section, which is 
similar to the early Universe thermally averaged s-wave dominated cross section {(J se ifv) = a. 
For s-wave annihilation with Xf ~ 20, and assuming that the candidate accounts for the 
whole of the DM, this is approximately {a se ifv) ~ 1.8 x 10~ 9 GeV _ . Note that (a se ifv) in 
the early Universe is equal to half of the annihilation cross section needed for non-identical 
particles with no asymmetry (A = 0) to make up all of the DM. Numerically, we find this 
cross section to be (c xx v) ~ 3.7 x 10~ 9 GeV -2 . 

For asymmetric DM, the annihilation rate is given by T asym = (a xx v)(p x p x /m x ), and 
we have assumed all along that p x + p x = pdm and n x = p x /m x , thus 



asym 



(°XX V ) 



Pdm Y x Y x 



(6.2) 



K (Y x + Y x y 

The ratio of the asymmetric DM annihilation rate to the standard symmetric annihilation 
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Figure 8. Same as in Fig. 6 but for scalar tensor models and {<J xx v) equal to 1 x 10 GeV - , 
8 x 10~ 6 GeV~ 2 , and 1 x 10~ 4 GeV~ 2 , respectively. The effect or the re-annihilation phase is clearly 
seen. 
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Figure 9. Comparison of numerical (blue) solutions for Y x (dashed lines) and Y x (solid lines) as a 
function of a; and analytical (red) solutions for Y x (x — > oo) and Y x (x ~ > oo) in the standard cosmology. 
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Figure 10. Same as Fig. 9 but for lO.a (left) kination and (<r xx v) = 1 x l(T 6 GeV _2 and lO.b (left) 
scalar-tensor cosmology and (cr xx v) = 8 x 10~ 6 GeV -2 . In both cases T tr =5 MeV. 
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Figure 11. 11. a (left) Contours of asymmetry A for which \ ano - X make up all the DM in the 
standard cosmology (red) , and in kination (blue) and scalar-tensor (green) models with T tr =5 MeV, 
for m x =100 GeV. 11. b (right) Evolution of Y x (dashed lines) and Y x (solid lines) for (<J xx v) = 5 x 
10~ 6 GeV -2 (indicated by the black vertical line in 11. a) in each of the three cosmological models. For 
this cross section A — 4.05 xl0~ 12 for the standard cosmology and kination models and A = 3.8 xl0~ 12 
for the scalar-tensor model. The evolution of Yeq for symmetric DM (A =0) is shown for comparison. 



In the standard cosmology, this ratio never exceeds one, but in the kination or scalar- 
tensor pre-BBN cosmological models, it may exceed one for a range of annihilation cross 
sections, as shown in Figs. 11, 12 and 13. This range of cross sections depends on the 
transition temperature Tt r at which the cosmology becomes standard. 

In Figs. 11, 12 and 13 we provide examples in which the asymmetric DM annihilation 
rate in a non-standard pre-BBN cosmological model is greater than the standard symmetric 
annihilation rate for m x = 100 GeV. We have chosen particular annihilation cross sections 
and transition temperatures in each case. In Fig. 11 the transition temperature is low, 



T tr =5 MeV and for (a 



xxv) = 5 x 10~ 6 GeV~ 2 , i.e. {&xx v ) I \ a seif v ) ~ IO 3 , for the scalar- 
tensor model we get Y x Yx/(Y x + Y X ) 2 ~ 10 -1 , so the annihilation rate is ~ IO 2 times greater 
than the standard symmetric DM annihilation rate. For the kination and standard models 
the present annihilation rate is negligible in this example, ~ 10~ 7 or smaller respectively 
than the standard symmetric rate. For kination, Y X Y X /(Y X + Y x ) 2 ~ 10~ 10 , and this ratio 
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Figure 12. Same as in Fig. 11 but with T tr — 40 MeV in the scalar-tensor model and (cr xx v) = 
9 x 1CT 7 GcV~ 2 for which A = 4.05 x 10~ 12 for the standard cosmology and A — 3.9 x 10~ 12 for the 
kination and scalar-tensor models. 
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Figure 13. Same as in Fig. 11 but with a high transition temperature T tr =3.35 GeV in both the 
kination and scalar-tensor models and {o x ^v) — 8.5 x 10~ 8 GeV~ 2 for which A — 3.95 x 10~ 12 for the 
standard cosmology and kination model, and A = 3.75 x 10~ 12 in the scalar-tensor model. 



is even smaller in the standard cosmology. Thus, in this example, the scalar-tensor model is 
easily distinguishable from the other two. 

For smaller annihilation cross sections we can select kination and scalar-tensor models 
which would produce almost identical effects in the annihilation rate, as demonstrated in Fig. 
12. In this figure, the transition temperatures, Tt r =5 MeV for kination and Tj r = 40 MeV 
for the scalar-tensor model, are chosen so that the contours of A for which \ an d X make up 



all the DM in both models are very similar. For (a 



xx l 



9 x 1CT 7 GeV -2 , both models give 
a present asymmetric DM annihilation rate ~ 10 times larger than the standard symmetric 
DM annihilation rate. If such a rate were detected, we could see that it is non-standard, but 
could not distinguish between the kination and scalar-tensor models chosen here. 

For asymmetric DM annihilation cross sections very close to the standard symmetric 
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DM annihilation cross section it is difficult to distinguish between any of the three pre-BBN 
cosmological models. In order to prevent the Universe from being overdense, the transition 
temperature of the non-standard cosmological models must be set very close to the standard 
cosmology freeze-out temperature. This means that % and x only evolve in the non-standard 
cosmological phase for a very short time after freeze-out. In Fig. 13, {a xx v) = 8.5 x 
10~ 8 GeV~ and the transition temperature is high, Tj r = 3.35 GeV for both kination and 
scalar-tensor models. All three models produce similar x relic densities, and the present 
asymmetric DM annihilation rate is ~ 10 times smaller than the standard symmetric DM 
annihilation rate for all of them. Thus the annihilation rate could not help us distinguish 
between the different pre-BBN cosmological models. 

These examples show that if the mass, approximate particle-antiparticle asymmetry, and 
annihilation cross section of a DM candidate is known, then if the present annihilation rate 
is detectable, one could possibly distinguish between these different pre-BBN cosmological 
models. A higher annihilation rate than predicted by the model of symmetric DM in the 
standard cosmology would imply the existence of a non-standard cosmology during pre-BBN 
epoch. 

6.1 Fermi Space Telescope bounds on non-standard asymmetric DM models 

In non-standard pre-BBN cosmological models the annihilation rate of asymmetric DM can 
be so large as to already be rejected by present bounds. The Fermi-LAT collaboration has 
published upper bounds [21] on the annihilation cross section of a symmetric DM candidate 
(a se ifv) < (o"v)Fcrmi as a function of its mass, from a combined analysis of Milky Way 
satellites, assuming that it constitutes all of the DM and annihilates only into particular 
final states, which correspond to limits on the annihilation rate in Eq. (6.1), rf±^P < rpermi- 
Here we use the XX ~^ ^fi an d XX ~^ bb modes, which impose the limit of (crv)Fermi equal 
to 8.80 x 10~ 25 cm 3 /s= 7.33 x 10" 8 GeV" 2 and 5.99 x 10" 26 cm 3 /s= 4.99 x 10~ 9 GeV" 2 to 
95% C.L. for m x = 100 GeV, respectively. Thus, under the same assumptions used by the 
Fermi-LAT collaboration, for our models, the ratio 



Lasym _ \ a xx v ) 2Y X Y X 



(6.4) 



1 Fermi \<TV /Fermi \*x ' x) 

must be less than 1. 

Fig. 14 shows (in pink) the regions of the DM asymmetry A versus asymmetric DM 
XX annihilation cross section in kination (Fig. 14. a) and scalar-tensor (Fig. 14. b) pre-BBN 
cosmological models already rejected by Fermi-LAT bounds [21] on the annihilation modes 
XX ~^ fJ>P> and XX ~~ ^ bb for m x =100 GeV, under the same assumptions made by the 
Fermi-LAT collaboration. These regions were found by varying the transition temperature 
of each of the two non-standard cosmological models between 5 MeV and 3.5 GeV. Each 
transition temperature provides a unique contour in the A versus cross section space that 
gives the right relic DM density. The Fermi-LAT bounds were then applied along each of 
these contours, dividing it into a section with acceptable annihilation rates and another 
with rejected annihilation rates. These contours were then combined to create the ruled out 
regions seen in Fig. 14. 

7 Summary and Conclusions 

In this paper we have studied asymmetric dark matter which decouples from the thermal bath 
during a non-standard pre-BBN cosmological phase. We have assumed that the dark matter 
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Figure 14. Regions of the DM asymmetry A versus asymmetric DM XX annihilation cross section in 
14. a (left) kination and 14. b (right) scalar-tensor pre-BBN cosmological models rejected by XX ~> MM 
(dark pink) and XX ~> bb (light pink) Fermi-LAT bounds [21] for m x =100 GeV, under the same 
assumptions made by the Fermi-LAT collaboration. 



particle asymmetry A is created prior to the period that we consider and that subsequently 
the only reactions which change particle and anti-particle numbers are pair annihilation into 
Standard Model particles and the inverse reactions; thus A is constant. We considered two 
non-standard pre-BBN models, kination and scalar-tensor models, in which the expansion 
rate of the Universe is respectively faster (~ T 3 ) and slower (~ T 12 ) than in the standard 
radiation dominated cosmology (~ T 2 ) until a transition temperature Xj r , which must be 
larger than about 4 MeV to preserve BBN and all the subsequent history of the Universe. 

We find that for a range of pair annihilation cross sections between about 5 x 1CP 6 
GeV -2 and the standard symmetric thermal dark matter cross section 1.8 x 10 -9 GeV -2 , the 
ratio of the relic abundance of the minority and majority dark matter components is highly 
dependent on the pre-BBN cosmology. For each annihilation cross section in this range, the 
smallness of this ratio in non-standard cosmological scenarios with particular Tj r values can 
be compensated and even overcompensated by the increased annihilation cross section so 
that the annihilation rate at present of asymmetric dark matter, contrary to general belief, 
could be even larger than that of symmetric dark matter in the standard cosmology. In 
the standard pre-BBN cosmology the annihilation rate of asymmetric dark matter is always 
negligible, because as soon as the annihilation cross section increases with respect to the value 
which would yield equal minority and majority relic densities, even by a factor of a few, the 
minority component density becomes exponentially smaller than the majority density. 

Thus, if the annihilation cross section of the asymmetric dark matter candidate is known, 
the annihilation rate at present, if detectable, could be used to test the Universe before Big 
Bang Nucleosynthesis, an epoch from which we do not yet have any data. 
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